A single-time two-point closure based on fluid particle displacements 
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A new single-time two-point closure is proposed, in which the equation for the two-point correla- 
tion between the displacement of a fluid particle and the velocity allows to estimate a Lagrangian 
timescale. This timescale is used to specify the nonlinear damping of triple correlations in the 
closure. A closed set of equations is obtained without ad-hoc constants. Taking advantage of the 
analogy between particle displacements and scalar fluctuations in isotropic turbulence subjected to 
a mean scalar gradient, the model is numerically integrated. Results for the energy spectrum are in 
agreement with classical scaling predictions. An estimate for the Kolmogorov constant is obtained. 



Two-point single-time closures are efficient and use- 
ful tools for studying homogeneous turbulence. Among 
existing single-time theories, the Eddy Damped Quasi- 
Normal Markovian theory (EDQNM), originally pro- 
posed by Orszag[l| (see also Leith [2[) for isotropic turbu- 
lence, has been used to investigate a broad range of fun- 
damental problems in turbulence. Examples are: scalar 
decay in isotropic turbulence ,3, 4], rotating turbulence 
stratified turbulence 0], magneto- hydrodynamics Q, rel- 



ative dispersion[8|, homogeneous shear [91 [1( 



isotropic 



turbulence with a mean scalar gradient [llL |l2j, decay 
of turbulence in a wall bounded domain [13| . compress- 
ible turbule nce fLUl , 2D turbulence [Hj], premixed flame 
propagation. [16| Without being exhaustive this list illus- 
trates the variety of problems that have been adressed 
using this single-time closure. EDQNM was also used 
to propose subgrid models for Large-Eddy Simulation 
fl7L fl8| . One common feature of all these works is 
that they rely on the original heuristic assumption of 
EDQNM: the presence of an eddy damping term. The 
eddy damping is an essential ingredient of the EDQNM 
closure, as it represents the nonlinear damping of the 
triple correlations, necessary to obtain inertial range 
spectra in agreement with classical Kolmogorov the- 
ory. It should be pointed out that in the case of weakly 
nonlinear wave-turbulence only the presence of the damp- 
ing is mandatory. [5] A convenient specification of its form 
is of primary importance as soon as strong turbulence is 
considered. In the case of EDQNM, the Eddy-Damping is 
specified by dimensional analysis and an ad-hoc constant 
is introduced in the model. For other closures like the 
Test Field Model proposed by Kraichnan, [2(| in which 
an auxiliary velocity field that does not respect the in- 
compressibility constraint is introduced, the damping is 
deduced from a more sophisticated analysis, but an ad- 
hoc constant still has to be introduced (of order unity in 
the case of the Test Field Model). 

More elaborate are the two-point two-time closures 
proposed by Kraichnan (DIA, for Direct Interaction Ap- 
proximation [2l|). obtained by perturbation techniques 
applied to the Navier Stokes equations. The Eulerian 
formulation of DIA, violating the principle of statisti- 
cal Galilean invariance and therefore being incompatible 



with a K 5 / 3 Kolmogorov inertial range, was reformu- 
lated in a Lagrangian framework. The Lagrangian His- 
tory DIA (LHDIA)Hi|, its abridged versions, [H HI] as 
well as the version of Kaneda [2J] are known to yield pre- 
dictions in agreement with Kolmogorov spectra without 
introducing any ad- hoc constant. The price one has to 
pay is the complexity of the models that depend on the 
entire Lagrangian history of the flow in contrast to the 
single-time theories. 

It is known that the EDQNM equations for isotropic 
turbulence can be obtained from two-time theories by as- 
suming an exponential decay of both the response func- 
tion and the two-time correlations. The eddy damping 
then corresponds to the inverse of the correlation time of 
the turbulent velocity field. In the framework of a La- 
grangian formulation of the theories, this correlation time 
has to be defined along fluid particle trajectories and in- 
deed a definition, corresponding to LHDIA for isotropic 
turbulence, can be found in Kraichnan [22I : 



r(K,t) 



1 E(K,t\s) 
o E(K,t) 



ds 



(1) 



with E(K,t) the energy spectrum and E(K,t\s) the La- 
grangian two-time energy spectrum (definitions are given 
below). 

In the EDQNM model, this time is modeled as a func- 
tion of E(K,t) and K, yielding by dimensional analysis: 



(r(K",i)) _1 = v(K,t) = ay/K 3 E(K,t) 



(2) 



Another variant, non-local in wave-number space, is the 
expression proposed in Pouquet et al. 



(r^i))- 1 = v (K,t) = X\ S 2 E{S,t)dS. (3) 



K 



At high Reynolds number, both formulations lead to 

T(K,t)~K- 2 / 3 (4) 

in the inertial range of the spectra, mimicking the scaling 
that the Lagrangian correlation timescale is expected to 



2 



follow in agreement with Kolmogorov theory. The con- 
stant a (or A) is specified to obtain the desired value of 
the Kolmogorov constant. 

The aim of the present letter is to derive a single-time 
closure that does not contain ad-hoc specification of the 
damping term nor adjustable constant. The eddy damp- 
ing timescale will be determined within the closure, the 
key element of the model being the use of equation |T|) . 

The Lagrangian two-time spectral tensor is defined by 

<£> tJ {K, t\s) = FT X _ X , [( Ui (x, t)uf(x', t\s))] (s < t), 

in which FT x _ x i denotes a Fourier transform with re- 
spect to x — x'. Uj(x',t\s) is defined as the velocity 
measured at time s within the fluid clement which passes 
through the point x' at time t. In isotropic turbulence 
the two-time energy spectrum E(K, t\s) is related to this 
tensor by the relation: 



E(K,t\s) 
AttK 2 



(6) 



with Pij (K) = (Sij -KiKj /K 2 ). E(K, t\s) is a key quan- 
tity as it appears in equation ([T|) giving the timescale that 
has to be specified in the closure. It is a two-time La- 
grangian quantity and is therefore difficult to evaluate in 
the framework of a one-time closure. However, only the 
time integral of E(K,t\s) is required to express r(K,t). 
This integral satisfies: 



AnK 2 



[ E{K,t\s)ds 
Jo 



Ui(x,t) / uf(x', t\s)ds 



(J) 



The integral of the Lagrangian velocity along the trajec- 
tory in ([7]) is the displacement of a fluid particle. Calling 
a,j this displacement: 



a,j(x,i) = Xj(x, t\t) — Xj(x, t\t = 0) 



uf(x, t\s)ds 



and 

Xj(x, t\t) = Xj 
([7]) can be written as: 
Pij{K) 



AttK 2 



E{K,t\s)ds = F Uia AK,t), 



(8) 



(9) 



(10) 



with 



F Uiaj (K,t) = FT SB - as ,[<u i (x,t)a j {x',t) >]. (11) 

In isotropic incompressible turbulence this quantity can 
be expressed as (Lumley [25j): 



F Uia AK,t) 



PjjjK) 
AttK 2 



F ua (K,t), 



(12) 



and therefore expression ([T]) can be recasted as a function 
of one-time Eulerian quantities only: 



r(K,t) 



F ua (K,t) 
E(K,t) 



(13) 



This expression forms the basis of the single-time two- 
point closure proposed in this paper. In the following, 
a way to obtain an expression for F ua (K, t) will be pro- 
posed that together with the equation for E(K, t) and 
relation (Ti"3"|) leads to a closed set of equations. 

From the Navier-Stokes equations and the equation for 
the displacement of a fluid particle: 



daj(x, t) 
dt 



Uj(x, t), 



(14) 



a one-time two-point closure for E(K,t) and F ua (K,t) 
in isotropic turbulence (and alternatively for &ij(K,t) 
and J- Ui a {K , t) for anisotropic turbulence) can straight- 
forwardly be written, by applying the Quasi-Normal ap- 
proximation and Markovian assumption and expressing 
the relaxation time of the triple correlations using (|13p . 
In the present letter, instead of deriving the evolution 
equation for F ua (K,t), we adopt a simpler approach, 
taking advantage of the analogy existing between the 
fluid particle displacement and an advected non-diffusive 
scalar field. This analogy will permit to express the 
closure model using only existing published equations. 
As pointed out by Batchelor, 26] a non-diffusive passive 
scalar in isotropic turbulence with a mean scalar gradi- 
ent obeys the same equation as the displacement of a 
fluid particle. Considering an isotropic turbulence ini- 
tially free from passive scalar fluctuations on which, at 
t = 0, a mean scalar gradient is imposed in an arbitrary 
direction, d<d/dxj, the interaction of the velocity field 
with the scalar gradient produces a scalar fluctuation 9 
governed by: 



ddjx, t) 
dt 



99 . , 



(15) 



Integrating (TT5|) over the Lagrangian trajectory of the 
fluid particle that arrives at time t at position x yields: 

59 

9(x,t) = - — (Xj(x,t\t) - Xj(x,t\t = 0)) 



—a^t) (16) 



and the correlation between 8(x',t) and Ui(x,t) can be 
expressed as: 

99 

< Ui(x,t)9(x' ,t) >= -- — < Ui(x,t)aj(x ,t) > . (17) 

OXj 

For isotropic turbulence without loss of generality, the 
direction of the gradient can arbitrarily be specified: for 
example x%. Equation (|17[) then leads to: 



<99 

F U3 g(K,t) = -—F U3m (K,t) 

Ox 3 



(18) 



3 



with F U3 g(K,t) defined as the scalar flux spectrum (see 
for example O'Gorman and Pullin [27|). or introducing 



Fua (K, t) as in equation (fl~3 



F U3 e(K,t) 



2 59 

3 9x 3 



The difference with the EDQNM model is that the 
eddy damping in the equations is not anymore heuristi- 
cally specified, but is calculated using relation (|20| . 
The equation for F U3 g(K,t) is: 



F ua {K,t). 



(19) 



F„ 



since isotropy implies that F uiai = F U2a , 

lp 
3 a ' 

One can also arbitrarily specify the value of the gra- 
dient, since the scalar equation is linear. Choosing 
dO/dx^ = —3/2 simplifies the formulation. With this 
particular value of the gradient, equation (|19[) simply ex- 
presses the identity between the spectrum of the scalar 
flux of the non-diffusive scalar and the spectrum of the 
velocity-displacement correlation. Hence it is straight- 
forward to use the EDQNM model proposed by Herr 
et al. [llj or Bos et al. for the scalar flux spectrum to 
calculate F U3 e(K, t) = F ua (K, t) (in this work we use the 
formulation of Bos et al.). The equations for E(K, t) and 
F ua (K,t) will be solved simultaneously to calculate the 
energy spectrum in isotropic turbulence and to evaluate 
the damping term using (|18[) that now takes the form: 



d_ 

dt 



vK 1 



F U3e {K, t) = P(K, t) + T» L e {K, t) + IL(K, t](25) 



which is the equation of reference (12) in the particu- 
lar case of a non-diffusive scalar. In (f25|) . P(K,t) is a 
term that in this case can be interpreted as the produc- 
tion of scalar flux by the mean scalar gradient, such that 
P(K,t) = -%§§-E(K,t) = E(K,t). The expressions for 
the nonlinear terms T^Jj(K, t) and H(K, t) are not repro- 
duced here (Equations (14) and (15) of reference (12)). 
These closed terms are exactly the same as in the case of 
the EDQNM model, except that the eddy damping is de- 
termined by equation (|20|) . More specifically, in relation 
(16) of reference (12), 



t^KPQ 



vt{K) + (i'(P) + n"(Q) + v{K 2 + P 2 ), (26) 



V (K,t) = (r(K,t))- 1 = 



E(K,t) 

FuAK,ty 



(20) 



// is replaced by rj and fj," is still zero as in Bos jl2i ]. 

The model is applied to the decay of isotropic turbu- 
lence by numerically integrating equations (|21[) and 
The energy spectrum is initialized by: 



The evolution equation for the energy spectrum is: 



E{K,0) = BK i e- 2K2 ' K i 



(27) 



d_ 

dt 



2vK 



E(K,t)=T NL (K,t), 



(21) 



T NL {K,t) 



with Kq = 1, Kl = 10 and B determined so that the ini- 
tial kinetic energy is equal to 1. The energy spectrum was 
evaluated during the period of self-similar decay. Spec- 
tra are shown at R\ = 150, 500 and 1500. The results 
in figure Q] show that a K~ 5 / 3 inertial range is obtained 
for the energy spectrum. The value of the Kolmogorov 
constant is estimated to be 1.73 as can be seen when the 
Pi Q) \ X V + z ) [K PE(P. t)E(Q, t) gpectrum is shown in compensated form. It has to be re- 
minded that in the case of the EDQNM closure this value 



in which the expression for the nonlinear transfer Tjvx 
is the classical single-time two-point closure expression 
(Orszag Q): 



-P E(Q,t)E(K,t) 



dPdQ i 
PQ 



where A is a band in P, Q-space so that the three wave 
vectors K,P,Q form a triangle. x,y,z are the cosines 
of the angles opposite to the sides K, P, Q of the tri- 
angle formed from K,P,Q. The characteristic time 
Q(K,P,Q) is defined as: 



Q(K,P,Q) 



1 - exp(-[i K pQ x t) 

f-KPQ 



(23) 



with 



VKPQ = ^(if 2 +P 2 +g 2 )+r ? (^,t)+7 ? (P,<)+ ?? (Q,<) (24) 



(2i&)not a prediction of the model but has to be specified by 
choosing the constant A. The value of Ck = 1.73 can be 
obtained with EDQNM by choosing the value A = 0.49 
in expression ([3]). A detailed comparison between the 
EDQNM model and the present closure deserves further 
attention. 

The results in figure Q] suggest that the present model 
yields a reasonable estimate of the Lagrangian timescale 
in isotropic turbulence. It would be useful to compare 
the present results to estimations of the timescale pro- 
vided by Direct Numerical Simulation using the method 
proposed in Lee et (see also Gotoh et al.^^) or to 

higher Reynolds number data provided by Large Eddy 
Simulation. [30l| 
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FIG. 1: Top: energy spectrum at R x = 150, 500 and 1500; 
bottom: compensated form. 
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